Absorption features caused by oscillations of electrons on the surface of a quark star 



(N 

o 

(N 
S3 

a 



X 

Or 
i 

O 
c3 



> 

SO 

(N 



X 



R. X. Xu\ S. I. Bastrukov 1 - 2 , F. Weber 3 , J. W. Yu 1 and I. V. Molodtsova 2 
1 State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing 100871, China 
2 Joint Institute for Nuclear Research, 141980 Dubna, Russia 

3 Department of Physics, San Diego State University, San Diego, California 92182, USA 

(Dated: January 10, 2012) 

If quark stars exist, they may be enveloped in thin electron layers (electron seas), which uniformly 
surround the entire star. These layers will be affected by the magnetic fields of quark stars in such 
a way that the electron seas would transmit hydromagnetic cyclotron waves, as studied in this 
paper. Particular attention is devoted to vortex hydrodynamical oscillations of the electron sea. 
The frequency spectrum of these oscillations is derived in analytic form. If the thermal X-ray 
spectra of quark stars are modulated by vortex hydrodynamical vibrations, the thermal spectra 
of compact stars, foremost central compact objects (CCOs) and X-ray dim isolated neutron stars 
(XDINSs), could be used to verify the existence of these vibrational modes observationally. The 
central compact object IE 1207.4-5209 appears particularly interesting in this context, since its 
absorption features at 0.7 keV and 1.4 keV can be comfortably explained in the framework of the 
hydro-cyclotron oscillation model. 

PACS numbers: 26.60.-c, 97.60.Jd, 97.60. Gb 



I. INTRODUCTION 

The spectral features of thermal X-ray emission are es- 
sential for us to understand the real nature of pulsar-like 
compact stars. Calculations show that atomic spectral 
lines form in the atmospheres of neutron stars. From 
the detection and identification of atomic lines in ther- 
mal X-ray spectra one can infer neutron star masses (M) 
and radii (R), since the redshift and broadening of the 
spectral lines depend on M/R and M/R 2 , respectively. 
Atomic features are expected to be detectable with the 
spectrographs on board of Chandra and XMM-Newton. 
No atomic features have yet been discovered with cer- 
tainty, however. This may have it origin in the very 
strong magnetic fields carried by neutron stars. An alter- 
native explanation could be that the underlying compact 
star is not a neutron star but a bare strange (quark mat- 
ter) star p], Q • The surface of such an object does not 
consist of atomic nuclei/ions, as it is the case for a neu- 
tron star, but of a sea of electrons which envelopes the 
quark matter. 

Strange stars are quark stars made of absolutely stable 
strange quark matter They consist of essentially 

equal numbers of up, down and strange quarks as well 
as of electrons |9i-flll]. The latter are needed to neutral- 
ize the electric charges of the quarks, rendering the in- 
terior of strange stars electrically neutral. Quark matter 
is bound by the strong interaction, while electrons are 
bound to quark matter by the electromagnetic interac- 
tion. Since the latter is long-range, some of the electrons 
in the surface region of a quark star reside outside of 
the quark matter boundary, leading to a quark matter 
core which is surrounded by a fairly thin (thousands of 
femtometer thick) sea of electrons [H, [t| [l2j . Due to the 
enormous advances in X-ray astronomy, more and more 
so-called dead pulsars are discovered, whose thermal ra- 
diation dominates over a very weak or negligibly small 



magnetospheric activity. The best absorption features 
(at ~ 0.7 keV and ~ 1.4 keV) were detected for the cen- 
tral compact object (CCO) IE 1207.4-5209 in the center 
of supernova remnant PKS 1209-51/52 (see Table [TJ) . Ini- 
tially, these features were thought to be associated with 



TABLE I. Dead pulsars (CCOs and XDINSs) with observed 
spectral absorption lines [13l. fTsL l23l] . with P the spin period, 
B the magnetic field (-Bio = -B/10 10 G) derived by magneto- 
dipole braking, T the effective thermal temperature detected 
at infinity, and Ea, the absorption energy. We do not list the 
B-fields of XDINSs for which the propeller braking could be 
significant because of their long periods. 



Source 


P/s 


-Bio 


fcT/keV 


£ a /keV 


RX J0822.0-4300 


0.112 


< 98 


0.4 




IE 1207.4-5209 


0.424 


< 33 


0.22 


0.7, 1.4 


CXOU J185238.6+004020 


0.105 


3.1 


0.3 




RX J0720.4-3125 


8.39 




0.085 


0.27 


RX J0806.4-4123 


11.37 




0.096 


0.46 


RX J0420.0-5022 


3.45 




0.045 


0.33 


RX J1308.6+2127 


10.31 




0.086 


0.3 


RX J1605.3+3249 






0.096 


0.45 


RX J2143.0+0654 


9.43 




0.104 


0.7 



the atomic transitions of ionized helium in a stellar at- 
mosphere where a strong magnetic field is present [3]. 
Soon thereafter, however, it was noted that these lines 
arc of electron-cyclotron origin [l5|. The spectrum of IE 
1207.4-5209 shows two more features that may be caused 
by resonant cyclotron absorption, one at ~ 2.1 keV and 
another, but of lower significance, at ~ 2.8 keV [Hj]. 
These features vary in phase with the star's rotation. 
Although the detailed mechanism which causes the ab- 
sorption features is still a matter of debate, timing obser- 
vations predict a rather weak magnetic field for this CCO, 



2 



in agreement to what is obtained under the assumption 
that the lowest-energy line at 0.7 keV is the electron- 
cyclotron fundamental, favoring the electron-cyclotron 
interpretation [13, Ell- Besides IE 1207.4-5209, broad 
absorption lines have also been discovered in other dead 
pulsars (listed in Table H} , especially in so-called X-ray 
dim isolated neutron stars (XDINSs), between about 0.3 
and 0.7 keV 0. 

In this paper, we re-investigate the physics of these 
absorption features. The key assumption that we make 
here is that these features originate from the electron seas 
on quarks stars rather than from neutron stars, whose 
surface properties are radically different from those of 
strange stars [1,0, HI]. Of key importance is the magnetic 
field carried by a quark star, which critically affects the 
global properties (hydrodynamic surface fluctuations) of 
the electron sea at the surface of the star. We study 
this problem in the framework of classical electrodynam- 
ics in terms of cyclotron resonances of electrons in weak 
magnetic fields, since the magnetic fields of dead pulsars 
are much lower than the critical field, B q = 4.414 x 10 13 
G, at which the quantization of the cyclotron orbits of 
electrons into Landau levels occurs. 



II. HYDRO-CYCLOTRON WAVES 
A. Governing equations 

For what follows we restrict ourselves to a discussion 
of the large-scale oscillations of an electron sea subjected 
to a stellar magnetic field. We will be applying the semi- 
classical approach of classical electron theory of metals 
and making use of standard equations of fluid-mechanics. 
The electrons are viewed as a viscous fluid of uniform 
density p = nm e (where n is the electron number den- 
sity) whose oscillations arc given in terms of the mean 
electron flow velocity 8v. This implies that the fluctua- 
tion current-carrying flow is described by the density of 
the convective current 5j = p e 8v, where p e = en is the 
electron charge density. The equations of motions of a 
viscous electron fluid are then given by [2(| 



The governing equation, Eq. ([IJ , can be represented as 



9Sw 1... . 
p— = - <5j x B + 77V 5v, 

at c 



(1) 



j = p e 8v, p = m e n, p e = en, (2) 

where e and n are the change and number densities of 
electrons, respectively. We emphasize that Sj stands 
for the convective current density and not for Ampere's 
j = (c/4-7r)V x JB, as it is the case for magneto- 
hydrodynamics. This means that the hydrodynamic os- 
cillations in question are of non-Alfven type. In Eq. ([TJ, r] 
denotes the effective viscosity of the electron fluid, which 
originates from collisions of electrons with the magnetic 
field lines at the stellar surface. It is worth noting that 
the cyclotron waves can be regarded as an analogue of 
the inertial waves in a rotating incompressible fluid, as 
presented in Eq. (III. 56) of Chandrasckhar's book [211 ] . 



+ w c [n B x dvj - 2/V ov = 0, 

eB Br) 
u c = , n B = — , v=-. 

m e c B p 



(3) 
(4) 



where ui c is the cyclotron frequency. In the Appendix, 
we show that the electron sea can transmit macroscopic 
perturbations in the form of rotational hydro-cyclotron 
waves which are characterized by the following dispersion 
relation, 



± uj c cos I 



1 



+ i- 



{vk 2 /uj) 



l-{vk 2 /uj) 2 \-(yk 2 juj) 



(5) 



where ui and k denote the frequency and wave vector of 
the perturbations, respectively. The Larmor radius of an 
electron in a strong magnetic field, tl — m e c 2 / (eB) oc 
B^ 1 , is very small for pulsar-like compact stars, and we 
neglect the viscosity term in the following analysis of 
the motion of collective electrons. In the collision-free 
regime, v = 0, the hydro-cyclotron electron wave is de- 
scribed as a transverse, circularly polarized wave whose 
dispersion relation and propagation speed are given by 



uj = ±uj c cos9 and V = ±(w c /fc) cos 9, 



(6) 



respectively. Here, 9 is the angle between the mag- 
netic field B and the wave vector k. If k || B one has 
uj = iuj c . In metals these kind of oscillations are ob- 
served as electron-cyclotron resonances. There are two 
possible resonance states, one for uj = uj c and the other 
for uj = —uj c . These resonances correspond to the two 
opposite orientations of circularly polarized electron cy- 
clotron waves. 



B. Hydro-cyclotron oscillations of electrons on 
bare strange quark stars 

We restrict our analysis to the collision-free regime of 
vortex hydro-cyclotron oscillations. Using spherical co- 
ordinates, equation ([3]) then takes the form 



^ V r r i 

— — +Lo c [n B x dv = 0. 
at 

Taking the curl of both sides of Eq. ([7]) , we obtain 
35 u> 



(7) 



dt 



uj c {n B ■ V) 5v, k = Vx <5v. (8) 



Let the magnetic field B be directed along the z-axis, so 
that in Cartesian coordinates ris = (0,0,1). We then 
have 



n r = cos 9, rig = — sin 9, = 0. 



(9) 



From a mathematical point of view, the problem can be 
considerably simplified if one expresses the velocity 5v, 
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which obeys the condition V • 8v — 0, in terms of Stokes' 
stream function, x(0,4>)- This leads to 



0, Svg 



1 dx(0,& 



rsittl 



dcf> 



The depth of the electron layer near the star is much 
smaller than the stellar radius so that r ~ R to a very 
good approximation. Equation ((HJ then simplifies to 



dSuj r 



-uj c 



n gSvg 



(11) 



with the radial component of the vortex given in terms 
of X, 



5uj r = 



dSvd 



V 2 = 



80 
1 



1 dSvg 

m9 d(j> 



= Vix(M),(12) 



89 



sin t 



1 



sin 2 90 2 



(13) 



Substituting Eq. ([12} and Eq. ([13]) into Eq. <(TTJ leads to 



'I 



<9x n 



(14) 



The fact that free electrons undergo cyclotron oscillations 
in the planes perpendicular to the magnetic field suggests 
that the stream function x c & n be written in the following 
separable form, 



X (0,0)=^(0)cos(<£±a;i). 



(15) 



The "+" sign allows for cyclotron oscillations which 
arc induced by the clockwise polarized wave, and the 
"— " sign allows for cyclotron oscillations induced by the 
count-clockwise polarized wave. Substituting (|15p into 
(fUJl leads to 



Vi MO) ± — MO) = 0. 



(16) 



In the reference frame where the polar axis is fixed, Eq. 
(Unj) is identical to the Legendre equation for the surface 
spherical function, 



VlP i (0)+£(£+l)P i (6) = O, 



(17) 



where Pe(cos9) denotes the Legendre polynomial of de- 
gree I. Hence, setting ip{9) = Pe(9) we obtain 



u±(i) = ± 



£(£ + !)■ 



eB 



£>1. 



(18) 



From this relation we can read off the frequency of a 
surface hydro-cyclotron oscillation of a given order I. 



C. Characteristic features of cyclotron frequencies 

Let us consider the spectrum of the positive branch 
u(£)=w+(£) of Eq. (HHJ), 



0J r 



1)' 



> 1. 



(19) 



TABLE II. Comparison between single-particle (Landau 
level) and hydro-wave results for the cyclotron frequencies 
(B12 = B/W 12 G, loi = lo(£ — 1), bj c denotes the cyclotron 
frequency). 



u(£ = 


1) u(l = 


2) u(l = 


3) fiwi/keV 


Landau level lj c 


2cj c 


3uj c 


II.6B12 


Hydro- wave oj c /2 


cj c /6 


oJc/12 


5.8B12 



From 



uj(£ + 1) 



> 1, 



(20) 



it follows that this ratio becomes a constant for £ > 1. 
Such a spectral feature is notably different from the one 
of electron-cyclotron resonances of transitions between 
different Landau levels. 

From the energy eigenvalues, E n , of an electron in a 
strong magnetic field, which are found by solving the 
Dirac equation (see Ref. [22|). one may approximate the 
value of E n for a relatively weak magnetic field, B <C B q , 
by 



mc + nhw c . n > 0. 



(21) 



Therefore, in the framework of a single-particle approx- 
imation, the emission/absorption frequencies, which are 
given by tkj(£) = E n+ ( — E n , should occur at 



lo{£) = £lo c , £>1. 



(22) 



Table ITU compares the results of Eq. (|22|) with the results 
of Eq. (|19|) obtained for the hydro-cyclotron wave model. 
Most notably, it follows that for the hydro-cyclotron wave 
model one obtains oj(£ = 2)/lj(£ = 3) = 2, in contrast to 
the cyclotron resonance model for single electrons which 
predicts this ratio for u(£ = 2)/u>{£ = 1) = 2.. 



III. IE 1207.4-5209 AND OTHER COMPACT 
OBJECTS 

As already mentioned in the Introduction, IE 1207.4- 
5209 (or J1210-5226) in PKS 1209-51/52 is one of the cen- 
tral compact objects in supernova remnants 23[, where 
broad absorption lines, near (0.7, 1.4) keV [14|. and pos- 
sibly near (2.1, 2.8) keV [n| were detected for the first 
time. The interpretation of the absorption feature at 
~ 2.8 keV is currently a matter of debate, in contrast to 
the feature at ~ 2.1 keV which is essentially unexplained. 
Intriguingly, an absorption feature with the same energy, 
2.1 keV, has also been detected in the accretion-driven 
X-ray pulsar 4U 1538-52 

For what follows, we assume that IE 1207.4-5209 is a 
strange quark star and that (some of) these absorption 
features are produced by the hydro-cyclotron oscillations 
of the electron sea at the surface of such an object. As- 
suming a magnetic surface field of B ~ 7 x 10 11 G and 
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TABLE III. The frequencies, ui(£), at which hydro-cyclotron 
oscillations occur for IE 1207.4-5209 with effective tempera- 
ture T — 0.2 keV, assuming a magnetic field of B ~ 7 x 10 11 
G. 

£ 1 2 3 4 5 6 

cj(£)/keV 4.2 1.4 0.7 0.4 0.3 0.2 



thus u>(£ = 3) = 0.7 keV, we obtain the oscillation fre- 
quencies shown in Table HTT1 A magnetic field of ~ 7X10 11 
G is compatible with the magnetic fields inferred for IE 
1207.4-5209 from timing solutions [lj] (9.9 x 10 10 G or 
2.4 x 10 11 G), since IE 1207.4-5209 shows no magneto- 
spheric activity and the P- value would be overestimated 
if one applies the spin-down power of magnetic-dipole ra- 
diation [25|, |26| . We note that the absorption feature at 
u>(i = 1) = 4.2 keV shown in Table HTT1 may not be de- 
tectable since the stellar temperature is only ^0.2 keV 
(see Table [J) , which will suppress any thermal feature in 
that energy range. 

Aside from IE 1207.4-5209, one may ask what would 
be the magnetic fields of other dead pulsars (e.g., radio- 
quiet compact objects) if their spectral absorption fea- 
tures would also be of hydro-cyclotron origin? Intrigu- 
ingly, the hydro-cyclotron wave model predicts magnetic 
fields that are twice as large as those derived from the 
electron cyclotron model if the absorption feature is at 
u>(i = 1) = lu c /2\ these fields could be ~ 10 times 
greater (see Table |TTJ) if the absorption feature is at 
u(£ = 2) = w c /6 or oj(£ = 3) = w c /12. The absorp- 
tion lines at (0.3 ~ 0.7) keV may indicate that the fields 
of XDINSs are on the order of - 10 10 to 10 11 G, if oscil- 
lation modes with I > 4 are not significant. 

As noted in [l5j |. unique absorption features on com- 
pact stars are only detectable with Chandra and XMM- 
Newton if the stellar magnetic fields are relatively weak 
(10 10 G to 10 11 G), since the stellar temperatures are only 
a few 0.1 keV. The fields of many pulsar-like objects are 
generally greater than this value, with the exception of 
old millisecond pulsars whose fields are on the order of 
10 8 G. Central compact objects, on the other hand, seem 
to have sufficiently weak magnetic fields (see Table HJ so 
that absorption features originating from their surfaces 
should be detectable by Chandra and XMM-Newton. Ar- 
guments favoring the interpretation of compact central 
objects as strange quark stars have been put forward in 
[27j , where it was shown that the magnetic field observed 
for some CCOs could be generated by small amounts of 
differential rotation between the quark matter core and 
the electron sea. 

Besides dead pulsars, anomalous X-ray pulsars (AXPs) 
and soft gamma-ray repeaters (SGRs) are enigmatic ob- 
jects which have become hot topics of modern astro- 
physics. Whether they are magnetars/quark stars is an 
open question (28[. In case that AXPs/SGRSs should be 
bare strange stars, the absorption lines detected from 
SGR 1806-20 could be understood in the framework 



TABLE IV. Comparison between the frequencies, to(£), and 
the absorption frequencies detected, c*J bs, for SGR 1806-20 
(with an assumption of B ~ 1.86 x 10 13 G). The data of w bs 
are from [2^]. Both cj(£) and aj b s are in keV. 



£ 


1 


2 


3 


4 


5 


6 7 8 


oj{£) 


108 


36 


18 


10.8 


7.2 


5.1 3.8 3.0 


Mobs 






17.5 ± 


.5 11.2 ± .4 


7.5 ±.3 


5.0 ±.2 - - 



of the hydro-cyclotron oscillation model. Assuming a 
normal magnetic field of B ~ 1.86 x 10 13 G so that 
fkj c /12 = 18 keV, one sees that the oscillation model 
predicts hydro-cyclotron frequencies which coincide with 
the observed listed in Table |TVj 

IV. SUMMARY 

In this paper, we study the global motion of the elec- 
tron seas on the surfaces of hypothetical strange quark 
stars. It is found that such electron seas may undergo 
hydro-cyclotron oscillations whose frequencies are given 
by u>(£) = uj c /[£(£ + 1)], where £ > 1 and w c the cy- 
clotron frequency. We propose that some of the absorp- 
tion features detected in the thermal X-ray spectra of 
dead (e.g., radio silent) compact objects may have their 
origin in excitations of these hydro-cyclotron oscillations 
of the electron sea, provided these stellar objects are in- 
terpreted as strange quark stars. The central compact 
object IE 1207.4-5209 appears particularly interesting. 
It shows an absorption feature at 0.7 keV which is not 
much stronger than the another absorption feature ob- 
served at 1.4 keV. This can be readily explained in the 
framework of the hydro-cyclotron oscillation model, since 
two lines with £ and £+1 could essentially have the same 
intensity. This is very different for the electron-cyclotron 
model, for which the oscillator strength of the first har- 
monic is much smaller than the oscillator strength of the 
fundamental. 

APPENDIX 

Here we derive the dispersion relation characterizing 
the propagation of hydro-cyclotron electron wave in the 
slab-geometry approximation. The governing equation of 
viscous electron fluid under the action of Lorentz force is 
given by 

d6v = El [6vxB]+rl y2 5 
at c 
which can be written as 

or 

+ u c [n B x 6v] - vV 2 5v = 0, (23) 

where 

eB B 77 

tu c = , n B = — , v=-, p = m e n, p e = en, 

m e c Bp 
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uj c is the cyclotron frequency, and rj stands for the effec- 
tive viscosity of an electron fluid originating from col- 
lisions between electrons. To make the problem ana- 
lytically tractable, we treat the electron sea as an in- 
compressible fluid and assume a uniform magnetic field. 
Equation (|23|) can then be written as 



(24) 



aSv „ 
—^r = -w c [n B x dvj + v V ov. 



Upon applying to Eq. (|24|) the operator V x , we arrive 
at 



at 



[V x Sv] = Lj c {n B -V)(5v-i/VxVxVx <5v(25) 



where V ■ Sv = 0, and 



at 



w c (n B -V)(5v-i/Vx Vxifu, (26) 



where Suj = V x 5~v. Considering a perturbation in the 
form of Sv = v' exp[i(kr — wt)], we have 



[k x <5v] = i— (n B ■ k)(5v + i— [k x Sv], (27) 
(l-i—\[kx6v]=i—(n B 'k)6v, (28) 

\ CJ / CO 



where (k ■ Sv) — 0. It is convenient to rewrite the last 
equation as 



[k x Sv] = i—(iiB ■ k) 5v 



I + i(vk 2 /uj) 



1 - {vk 2 /ujf 



(29) 



Multiplication of both sides of Eq. (j29|) with k leads to 

1 + i(vk 2 /uj)~ 



jSvk = iuj c (riB ■ k) 



I - {uk 2 /uj) 2 



[k x <5v].(30) 



Inserting the left-hand-side of Eq. (|29| into the right- 
hand-side of Eq. (|30p gives 



; ("B ■ k ) 

k 2 



2 r 



1 + i(i>k 2 /w) 



or 



(ng ■ k) 

fc 



f - (^fc 2 /^) 2 
I + i(yk 2 jbj) 



(31) 



1 - (i/fc 2 /^) 5 



which is Ec;. ([5])- 
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